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m _ 

| We use a homogeneous sample of about 1100 optical and radio rotation curves 

(RCs) and relative surface photometry to investigate the main mass structure proper- 
ties of spirals, over a range of 6 magnitudes and out to ^ 1.5 and 2 optical radii (for the 
optical and radio data, respectively) . We definitely confirm the strong dependence on 
| luminosity for both the profile and the amplitude of RCs claimed by Persic & Salucci 

(1991). Spiral RCs show the striking feature that a single global parameter, e.g. lu- 
minosity, dictates the rotation velocity at any radius for any object, so unveiling the 
existence of a Universal RC. At high luminosities, there is a slight discrepancy between 
the profiles of RCs and those predicted from the luminous matter (LM) distributions: 
this implies a small, yet detectable, amount of dark matter (DM). At low luminosities, 
the failure of the LM prediction is much more severe, and the DM is the only relevant 
mass component. We show that the Universal RC implies a number of scaling proper- 
ties between dark and luminous galactic structure parameters: (a) the DM/LM mass 
ratio scales inversely with luminosity; (b) the central halo density scales as L -0 7 ; (c) 
the halo core radius is comparable to the optical radius, but shrinks for low luminosi- 
ties; (d) the total halo mass scales as L 5 . Such scaling properties can be represented 
as a curve in the (luminosity)-(DM/LM mass ratio)-(DM core radius)-(DM central 
■ density) space, which provides a geometrical description of the tight coupling between 

the dark and the luminous matter in spiral galaxies. 



1 INTRODUCTION 

It is well known that rotation curves (hereafter RCs) of spiral galaxies do not show any Keplerian fall-off which implies, as 
a most natural explanation, the presence of an additional and invisible mass component (Rubin et al. 1980; Bosma 1981b; 
Faber & Gallagher 1979; see also Ashman 1992). In more detail, the profiles of the RCs imply that the distribution of light 
does not match the distribution of mass: in fact, in each galaxy the local M/L ratio increases with radius by up to a factor 
~ 10 3 . Moreover, decreasing the galaxy luminosity, the light is progressively more unable to trace the radial distribution of the 
dynamical mass (e.g.: Persic & Salucci 1988, 1990a,b, 1991; Broeils 1992b). This is evident in the Universal Rotation Curve of 
spirals, pioneered by Rubin and collaborators (e.g., Rubin et al. 1985; see also Burstein & Rubin 1985 for the related concept 
of 'mass types') and derived by Persic & Salucci (1991; hereafter PS91) from a sample of 58 late-type spiral RCs. PS91 found 
that 

V(R) ~ 200 (£\ [1 + f (L b / L* b ) km s- 1 0.5< R/R ovt < 1 (1) 



© 0000 RAS 



2 M.Persic, P.Salucci and F.Stel 



with a ~ 0.25, Rm = 2.2 Rd logLJ = 10.4 {Lb denotes B luminosities)^, and fiLs / 1 L* B ) a function linear in log(Ls j 'L B ) 
which takes the values of 0.5 and at log(Ls / L* B ) = —1.2 and 0.4, respectively. The main features of eq.(l) are: (a) the linear 
dependence on radius, and (b) the luminosity dependence of the slope 

v = %v\*~ = 9{Lb/l ' b) 

where g(L) is a decreasing function of luminosity ranging between —0.1 and 0.6 (see PS91). A similar luminosity dependence 
of the RC shapes is seen also at outer radii (i.e., for R op tlS R^, 2R opt ): the RCs of low-luminosity galaxies continue to rise and 
then flatten out asymptotically, while the RCs of high-luminosity objects show a distinctive drop just outside R op t, followed 
by an asymptotically flat regime with Vac < V(R op t) (Salucci & Frenk 1989; see also Casertano & van Gorkom 1991). 

This dependence of the RC shapes has a pivotal importance in the dark matter (DM) issue because it implies that the 
structural properties of dark and luminous matter are connected (see Appendix A), probably as a consequence of the process 
of galaxy formation itself. In addition, the study of RC profiles can provide a straightforward test for power-law perturbation 
spectra (Moore 1994), and so it is a crucial aspect in formulating theories of galaxy formation (e.g., Evrard et al. 1994; Navarro 
& White 1994; Flores et al. 1993; Navarro, Frenk & White 1996). 

On the other hand, over the last few years the amount of available data has increased by more than one order of magnitude 
(e.g., Broeils 1992b; Amram et al. 1992, 1994; Schommer et al. 1993; Persic & Salucci 1995, hereafter PS95). For instance, 
the PS95 sample (967 curves) is a factor ~15 bigger than that used in PS91. It is worthwhile, then, to make a definitive 
assessment of our previous results (e.g., Persic & Salucci 1988, 1990b, 1991) and of their interpretation in the light of the 
recent theoretical work. In detail, the aim of this paper is to derive the Universal Rotation Curve of spiral galaxies and to 
investigate the main properties of the dark matter distribution. 

The plan of this paper is the following: in Section 2 we describe the selection criteria and the procedure used to build the 
RC samples; in Section 3 we introduce two model-independent DM indicators and we study the DM systematics; in Section 
4 we derive the Universal Rotation Curve; in Section 5 we discuss some implications of our results for current scenarios of 
galaxy formation. Finally, in Section 6 we comment on our main results and draw conclusions. 



2 SELECTION CRITERIA AND RC SAMPLES 

In order to evaluate the ability of a RC to probe the distribution of matter out to the outskirts of the distribution of the 
luminous matter (LM), we set a suitable reference scale: R op t, the radius encompassing 83% of the total integrated light. For 
an exponential surface brightness distribution, I(r) oc e~ r/RD , we have R op t = 3.2 Rn-f\- Therefore, we can consider that the 
stellar disc has a size ~ R opt and that no appreciable amount of LM exists farther out. Obviously a RC not extended out to 
R op t gives little information on the underlying gravitational potential (see Lake & Feinswog 1989). Let us now introduce the 
samples of RCs: 

• Sample A. It includes 131 individual RCs with a reliable profile out to R op t- To ensure this, we have set the following 
selection criteria. In order to be included, an RC must: (a) extend out to > R op t; (b) have at least 30 velocity measurements 
distributed homogeneously with radius and between the two arms; and (c) show no global asymmetries or significant non- 
circular motions: the profiles of the approaching and receding arms must be compatible. For 21-cm RCs we additionally 
require that (d) the beam-size be < 1/3 R op t- Let us notice that criterion (b) is also required to apply criterion (c). An 
RC fulfilling these requirements is generally smooth, with a well-defined profile whose average slope can be estimated with 
sufficient precision. On the other hand, an RC which fails the above criteria has an ill-defined radial profile and does not 
faithfully represent the underlying gravitational potential. In Appendix D we report the relevant photometric and kinematical 
quantities (Table 1) alongside their references (Table 2). 

• Sample B. The RCs of PS95 are co-added to form 22 synthetic curves, arranged by 11 luminosity bins and 11 velocity 
bins. Since this procedure averages out most of the observational errors and non-axisymmetric disturbances present in the 
individual RCs, only criterion (c) is applied and this only to the most severe cases (67 RCs). Instead, we require that RCs be 
extended at least out to 0.8 R op t [i.e., a relaxation of criterion (a)], which leads to the elimination of 284 RCs. The final set 
comprises 616 curves. We derive the synthetic RCs of Sample B as follows. 

(i) We normalize each RC to Ves, the velocity at the radius encompassing 65% of the integrated light: for a self-gravitating 
exponential disc, Ves corresponds to the radius (2.2 Rd) where the circular velocity peaks (see Freeman 1970). With this 
normalization, at any radius it can immediately be seen whether the light traces the dynamical mass, by looking at the 
amplitude and at the profile of a given rotation curve. 

• Rd is the disc exponential length-scale; R op t = 3.2 i?o is the optical radius (see below), 
t A value of the Hubble costant of Ho = 75 km s — 1 Mpc -1 is used throughout this paper. 

T For a Freeman disc this corresponds to the de Vaucouleurs 25 S-mag/arcsec 2 photometric radius. 
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Figure 1. Synthetic rotation curves for Sample B arranged by luminosity. Galactocentric radii are normalized to R op t, the radius 
encompassing 83% of the total / luminosity. The last panel shows the rotation curve predicted for a pure self-gravitating exponential 
thin disc. 



(ii) We divide the 616 RCs into 11 luminosity intervals spanning the whole /-band luminosity range —16.3 < Mi < —23.4, 
and into 11 velocity-amplitude intervals spanning the range 40 < V a /(km s~ ) < 340 (V a is taken from Mathewson et al. 
1992). Each interval, of central luminosity logL (velocity amplitude V a ), typically contains 70 RCs. Within each interval, 
the normalized velocity data are co-added to obtain the corresponding raw synthetic RCs, V(R;\ogL) and V(R; V a )- Each 
of these has about ~ 1500 velocity measurements, and represents the (average) RC for a spiral of luminosity log L (velocity 
amplitude V a ). 

(in) Finally, we form the smoothed synthetic RCs by averaging the data in radial bins of size 0.1 Ropt • 

In Fig.l (and Fig. 14 in Appendix B) we show the luminosity (velocity) sequence of the synthetic RCs. We find that, in 
the region 0.4^ R/R op t^ 1.1, the RCs are essentially straight lines, justifying the procedure of Persic k, Salucci (1988, 1990b) 
in deriving their slopes. The present analysis, by co-adding several hundred RCs, is negligibly affected by: a) errors in the 
assumed galaxy distances and/or inclinations; b) non-axisymmetric disturbances in the RCs; and c) any kind of observational 
errors. Moreover, we do not assume here a specific functional form to represent the RC. As a result, we confirm (Persic & 
Salucci 1990b) that the inner slope is a strong (inverse) function of luminosity: steeply rising curves and slightly falling curves 
are found at the two opposite extremes of the luminosity range, V ~ 0.5 at Mj = — 17 and V ~ —0.1 at Mi — —24. 

The co-addition procedure averages over empirical properties of the galaxies such as Hubble type, strength of bulge, colour, 
spiral design, and the presence or absence of bars. Since these quantities show, at most, a very loose trend with luminosity 
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(de Vaucouleurs et al. 1991; Biviano et al. 1991; Elmegreen & Elmegreen 1990), we leave the study of the dependence of the 
Universal Rotation Curve on galaxy morphology to a separate paper (see, however, Section 5) . 

3 DM INDICATORS: SYSTEMATICS OF DM 

The existence of DM in spiral galaxies is evident from the above considerations. The next step, i.e. to derive the mass 
structure of the spirals, requires assumption of a mass model in order to analyse the rotation curves appropriately. However, 
Persic & Salucci (1988, 1990b) have shown that the main features of the dark matter distribution can be obtained more 
straightforwardly by comparing the light profile and the mass profile (i.e. the RC shape): the inability of the luminous 
matter to trace the distribution of gravitating matter becomes the signature of a dark component. In detail, given a sample 
with photometric and kinematical data, at R opt we compute: (a) the slope V; um = (d\ogVi um /dlogR)n opt , where Vi um (R) 
is determined by the self-gravity of the visible matter derived from the surface brightness profile [with no assumption on 
(M I L)i um ]; and (b) the actual slope of the rotation curve V = (d\ogV / dlogR) n opt . Vj um and V are computed using the 
average slopes of Vi um (R) and V(R) in the range (0.6 — 1) R op t (see, e.g., Persic & Salucci 1990b). 

Let us now introduce the quantity V — V; um which compares the distribution of gravitating matter with that of the 
luminous matter, ensuring a model-independent detection of a mass discrepancy inside R op t- In fact, the light traces the mass 
as long as V — Vi um < 0.04 (where the r.h.s. term represents the typical uncertainty in V), while a strong disagreement 
between the two slopes (i.e., V — Vi um lt 0.2) implies the presence of a non-luminous component. The difference has an 
immediate interpretation because, given the similarity of the surface-brightness profiles of normal spirals (see Fig. 12) and the 
relative smallness of the bulge and gas contributions, V; um does not depend on luminosity. In fact we have estimated (see 
Appendix A) that 

Vi„ m = -0.24 ±0.03. (2) 

Therefore the quantity V + 0.25 probes a mass discrepancy inside R op t that we will interpret in terms of a dark component. 
The values of V are reported in Table 3a for the 131 individual RCs of Sample A, and in Tables 3b,c for the 11 + 11 synthetic 
RCs constructed from Sample B. 

In our samples many rotation curves are extended beyond the optical size, out to ~ 2R op t. At these radii the disc and 
bulge contributions have virtually reached the Keplerian regime, V oc r~ 1//2 , and such fall-off must be apparent in the RC 
unless the stellar component is negligible. Let us introduce a mass-discrepancy indicator that tests the observed fractional 
variation in velocity, 8 = [V(2R opt ) — V(R op t)]/V(R op t), with the variation 8i um that we would expect if the stars and the 
gas actually traced the mass. Consequently, the quantity 8 — 8i um indicates the level of mass discrepancy inside 2R opt . The 
outer gradients 8 are mostly estimated from HI data: in Sample A we found a subset of 27 RCs extended to about 2 R opt (see 
Table 1 in Appendix D). In addition, about 200 optical RCs of Sample B have measurements farther out than R op t' these 
data, though relatively sparse, well supplement the HI data (see Appendix B). Then we can build the largest sample for which 
the haloes' RCs can be traced directly from the observed RCs, in that the latter are extended out to radii where the disc 
contribution is certainly negligible. The gradient 8i um can easily be estimated from surface photometry and the HI density 
distribution: we find that this quantity is uncorrelated with galaxy luminosity or velocity amplitude (see Appendix A) : 

Slum = -0.25 ±0.02. (3) 

The mass-discrepancy indicator is then simply given by 8+ \. Values of 8ft imply for R ~ 2R op t a strong decoupling between 
mass and light. 

With these indicators we now examine the presence of dark matter in spirals. 



3.1 Dark Matter inside R opt 

First, we quantify in a simple way the trend between RC profiles and luminosity. In Fig. 2 we plot V versus luminosity (velocity) 
for the objects of sample A. An inverse proportionality between slope and luminosity (velocity) emerges very strongly: the 
probabilities of a random occurrence are < 10~ 5 %. The constant of proportionality is —0.39 ± 0.02 ( — 1.26 ± 0.06), the linear 
correlation coefficient is 0.83 (0.87), and the scatter is <7v ; l b = 0.114 (cr V; v(R opt ) — 0.099). Remarkably, the co-added sample 
shows this very same trend but with an even higher level of statistical significance and less noise (see Fig 2). This suggests 
that part of the scatter around the above relationships arises from observational errors or kinematical disturbances, obviously 
suppressed by the the co-addition procedure: at a given luminosity (velocity amplitude), the intrinsic variance of V is then 
ctv~ 0.05, one order of magnitude smaller than the range of V spanned by spirals. From the combined two samples we obtain, 
with probability of random occurrence < 10~ 6 , 

V = 0.10 - 0.36 x log-^- (4a) 
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TABLE 3a. Sample A: Inner Gradients 
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-0.07 
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116-012 


0.49 
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-088 


0.14 
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0.41 
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0.02 
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0.29 


121-006 


0.53 
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■030 


0.05 
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0.28 


N 
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0.60 


N 7536 


0.08 


123-023 


0.41 


439- 


■018 


-0.12 
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0.39 


N 


3145 


0.04 
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0.03 


141-020 


0.00 


439- 


■020 


0.09 


N 598 


0.36 


N 


3198 


0.19 


N 7593 


0.17 


141-034 


-0.04 


444- 


■047 


0.45 


N 628 


0.00 


N 


3200 


0.00 


N 7606 


-0.06 


184-051 


-0.04 


444- 


■0S6 


0.19 


N 697 


-0.03 


N 


3223 


0.00 


N 7631 


0.19 


215-039 


0.25 


445- 
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0.05 


N 753 


0.02 


N 


3992 


0.02 


N 7793 


0.56 


235-016 


0.08 


446- 


■044 


0.47 
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-0.07 


N 


4013 


0.00 
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0.23 


240-011 


0.00 


481- 


■002 


0.30 
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0.00 


N 
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0.22 
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0.00 
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0.26 
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4236 


0.41 
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0.11 
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N 
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N 
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-014 


0.36 
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0.04 


N 
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0.00 


27-008 


0.27 
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0.20 
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-020 


0.52 
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0.52 


N 
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-0.04 
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-0.08 
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0.61 


566- 


■022 


0.26 


N 1832 
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N 


5585 


0.41 


40-012 


0.09 
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0.24 
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■009 


-0.07 


N 2336 


-0.08 


N 
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0.25 


41-009 


0.24 


350-023 


0.00 
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0.39 


N 2403 


0.33 


N 
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0.04 
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0.20 
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0.04 
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0.00 


N 
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N 
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0.26 
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0.32 


376-002 
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0.18 


N 
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-0.16 
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0.04 


379-006 


0.29 









TABLE 3b. Sample B: Inner Gradients (Luminosity Bins) TABLE 3c. Sample B: Inner Gradients (Velocity Bins) 

Mj Vio-y V a (km/s) V±cr v 



-18.50 





.53 


±0. 


.04 


76 





.55 ±0.04 


-19.44 





.48 


±0. 


.03 


103 





.50 ±0.03 


-19.97 





.41 


±0, 


,03 


120 





.40 ±0.03 


-20.50 





.35 


±0. 


.03 
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.34 ±0.02 


-20.88 





.30 


±0. 


.03 


154 





.27 ±0.02 


-21.25 





.25 


±0. 


.03 
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.22 ±0.02 


-21.57 





.20 


±0. 


.03 


187 





.18 ±0.02 


-21.96 





.15 


±0 


.02 


208 





10 ±0.02 


-22.25 





.10 


±0. 


.03 


229 





07 ±0.04 


-22.64 





.03 


±0. 


.02 


254 


-0. 


02 ±0.03 


-23.19 


-0 


07 


±0 


02 


300 


-0. 


11 ±0.04 



V = 0.10 - 1.35 x log-^jg, (46) 

where the uncertainties on the coefficients are less than 8%. We remark that the inner slope V turns out to be a tight and 
continuous relationship of luminosity, as found in PS91. In disagreement with the claim of Flores et al. (1993) (see Section 
5), the RCs are very steep at low luminosities (velocity amplitudes), and become gradually flat — and eventually slightly 
decreasing — at higher luminosities (velocities) . 

Recalling that V; um ~ —0.25, relationship (4) has two straightforward implications: (a) in every galaxy the LM is unable 
to account for the observed kinematics out to ~ R op t\ and (b) this inability is marginal (though clearly detectable) at high 
luminosities, and becomes progressively more serious with decreasing luminosity. 

We estimate the minimum mass, (1 — /3 max )G~ 1 V 2 (R op t)Ro P t, of the dark component residing inside R opt by assuming 
the maximum- disc hypothesis: 

Pma* = 5(1_ V ) ~ | + Jlogi- (5) 

(for details see Persic & Salucci (1990c). Then, in the biggest (brightest) spirals the dark matter comprises at least ~ 15% of 
the total mass, while in the smallest (faintest) ones this fraction increases to ~ 85%. 
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Figure 2. The RC slope in the region (0.6 — 1) R op t versus absolute magnitude (a) and versus (log of) rotation amplitude (b) 
[V op t = V(Ropt)]- The filled circles represent the individual curves of Sample A; the empty squares represent the synthetic curves of 
Sample B. The conversion between J and B magnitudes has been obtained by cross-correlating Sample B with the RC3 catalog (dc 
Vaucouleurs et al. 1991): Mb = —0.38 + 0.92M/ (129 objects; this relation will be used throughout the paper). 
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Figure 3. The fractional variation <5 of the rotation velocity between R op t and 2 R op t versus absolute magnitude (a) and versus the 
(log of) rotation velocity at R op t (b). Filled circles and empty squares denote, respectively, the individual objects in Sample A and the 
extended synthetic RC's from Sample B (given in Fig.B2 of Appendix B). 



3.2 Dark matter inside 2R op t 

In Fig. 3 we plot S versus luminosity and versus velocity for the very extended HI RCs of Sample A: a tight outer slope- 
luminosity (velocity) relationship emerges (the probability of a chance occurrence is less than 10 -4 ). Also shown in Fig.3 are 
the estimated 5s for the 5 synthetic extended RCs from Sample B (see Fig. 15 in Appendix B): the trend is clearly consistent 
with that of individual RCs from Sample A. All these data are very well fitted by: 

8 = -0.05 - 0.16 log-^-, (6) 

S = - 0.05 - 0.47 log^ , (7) 

(with uncertainties of 20% in zero point and 10% in slope), which confirms and strengthens the results of Salucci & Frenk 
(1989) and Casertano & van Gorkom (1991) by means of the consistent evidence of about 100 RCs. 

By recalling that 5i um ~ —0.25, we can state that, at ~ 2R opt , the luminous matter is totally unable to account for the 
observed kinematics. 
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Figure 4. The Universal Rotation Curve of spiral galaxies. Radii are in units of R op t- 



3.3 Implications 

The main implications of the above results can be summarized as follows. 

• The fraction of DM increases with radius: in fact, if (1 — (5) is the fraction of dark matter inside R op t, at 2 R opt this fraction 
has increased to ~ 1 — 2 (i+s)' 2 ■ 

• The maximum disc mass is M™£f ~ 2.6 x 10 11 /if 5 1 M Q , while the minimum disc mass is found to be M™m ~ 7 x lO 8 /if 5 1 M . 

• The fraction of dark matter increases, at any normalized radius R/R opt , with decreasing galaxy luminosity. 

• The dark component is much less concentrated than the visible one: in fact, since puM oc V 2 /R 2 and plm oc R~ 2 e~ 3 2 fl ' Ropt , 
the dark-to- visible density ratio at 2R opt is 15-30 times larger than the value at R opt . 

Thus the ensuing picture of a spiral galaxy consists of a stellar disc embedded in a dark component with a very different 
density profile. Spirals of high luminosity (Ms < —21) have a small dark matter content inside R op t - consequently, their 
RCs show the kinematical signature, 8i um < 8 < 0, of the transition from an inner, LM-dominated region, to an outer, DM- 
dominated region (Salucci & Frenk 1989). Conversely, low-luminosity (Mb^ — 20) objects are dominated nearly everywhere 
by the DM component, 8 > 0, so that the transition occurs well inside R op t where Vi um has still an increasing or flattish 
profile not dissimilar to that of the DM, and therefore no feature marking the transition between the two regimes appears in 
V(R). 



4 THE UNIVERSAL ROTATION CURVE 

The luminosity specifies the entire rotation field V(R, L) of spirals. From samples A and B we obtain, and plot in Fig. 4 
the synthetic rotation curve V(R;L) [and V(R; V a )\ as a function of radius and of galaxy luminosity (velocity). These are 
constructed by combining the co-added rotation curves with the slope- velocity-luminosity relationships (6), (7). The rms errors 
of the 11x2 co-added RCs are typically 2 ± 1% inside R op t, while beyond R op t they increase by up to a factor of 2. 

Crucially, V(R; L) has a very small cosmic variance. In fact, the percentage residuals A between the individual RCs, 
Vi(R; Li) (given in Table 2 of PS95), and the averaged RCs, V(R; L), are generally within the estimated observational errors 
(2 — 4 per cent). In Fig. 5 we show the comulative distribution of A: only about ~ 1/6 of points has a variance that cannot 
be simply accounted for by observational errors. However, also in these cases, the variance is rather small (~ 15%), and 
furthermore may arise from a neglected dependence of RCs on Hubble Type, or from non-gaussian observational errors, 
rather than indicating actual violations of the Universal Rotation Curve paradigm. Thus, as in PS91, we here claim the 
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Figure 5. The cumulative distribution of the relative differences between individual and averaged RCs (all luminosity bins are used). 

existence of a Universal Rotation Curve, V(R;L) (given in numerical form by 'Table URC |§] and shown in Fig. 4). (Sa-type 
galaxies will be investigated in a separate paper.) 

We now look for an analytical form for V(R, L), that we define as Vurc{R', L). First, we stress that a very simple function, 
i.e. two straight lines smoothly joining at R op t with luminosity-dependent slopes, describes well the typical profile of RCs (see 
PS91): the inaccuracies are J$ 5% (see Persic, Salucci & Stel 1996). Here we proceed a step further by fitting V(R;L) with 
a function that is still simple, but is nevertheless able to: (a) describe the innermost parts of the RCs; and (b) pinpoint the 
separate contributions of the dark and luminous components. In detail: 

(i) The contribution from the stellar disc can be written, for 0.04R op t < R < 2R opt , as (see Appendix C): 

Vi( X ) = V\R opt ) (3 (j2 + o -782)1 .43 , (8) 
where x = RIR op t and 8 = \Vd(R op t) IV (R op t)] 2 ■ For an exponential thin disc, B = ,A', 1 b GM - i L b — I 

* \ftopt) ft opt 

(ii) The contribution from a dark halo can be well represented by: 

Vtf(s) = V\R opt ) (1 - /J) (1 + a 2 ) (a / + 2 Q2) , (9) 

with a being the 'velocity core radius' (measured in units of R op t)- We normalize Vurc by setting VuRc(Ro P t) = V(R op t\ L). 
Therefore, (1 - B) = v ^f k J . 

Then we fit V(R; L) with the universal rotation curve Vurc given by: 

VuR C {R;B,a) = [vi(R) + Vt(R)] 1/2 . (10) 

We find that, with 

B = 0.72 + 0.44 log-^- (11a) 



a = 1.5 , (116) 

the universal rotation curve Vurc(R', L) reproduces the observed rotation curves V(R; L) to well within their rms errors (see 
Fig.6). 

The introduction of a bulge component should further improve the agreement between data and fit, by reducing the small 
discrepancies between eq.(10) and the data for very small radii in some synthetic RCs of high luminosity (see Fig.6). In Fig.6 
the RCs, V( R R - ; L), are gauged to the reference frame of the luminous matter: the differential effect of the presence of DM 
emerges clearly as an increasing failure of the LM to account for the observed rotation curves as luminosity decreases. 

We now set the RCs in the reference frame of the dark matter, whose distribution can be characterized by -R200, i.e. 

§ In /pub/psrot at 147.122.2.158 via anonymous FTP. 

In spherical symmetry. Otherwise, Md m represents the effective mass. 
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Figure 6. Best two-component fits to the universal rotation curve (dotted line: disc; dashed line: halo). The URC beyond R op t is built 
by linear extrapolation according to eq. (6). Notice that the extent of the RCs and the smallness of their rms errors limit the uncertainties 
on the parameters f3 and a to about 10% and 5%, respectively. 



the radius encompassing a mean halo overdensity of < 5p/p >= 200. To smooth the density field we use a top-hat filter: 
< p >r 20 o = Mdm (-R20o)/(§7r-Rioo); with the halo mass Mdm(R) = J 4irr 2 pn{r) dr. (In this formalism the 'central halo 



3 V ( 

density' is given by lim_R^o ^JqjP ano - does not depend on the local density, 
density of the Universe is p — 3Hq/(8itG) |J| then R200 is obtained by solving 



oc aa ' D M ^ dR .) Recalling that the mean mass 



/5ph\ 
\ P 1 1 



H 2 



V 2 



R 2 



200. 



^200 U L - 1 ^200 

with V(R) = V UR c(R) for 7? < 2R opU and V(R) = V URC (2R opt ) for R > 2R opt . 
(12), can be obtained either directly from (9) and (11) or by using the relation 

/ r \ - 4 - 

Vi(R 2 



The quantity V 2 



(12) 



V 2 , appearing in 



V 2 (R op t)R 20 



R. 



opt 



(13) 



(see Persic & Salucci 1990b): the two estimates are in very good mutual agreement. [In any case, let us remark that Mlm « 
M200, with M200 = Mdm{Rwo)\ In Fig. 7 we plot V{ R ^ i L): we realize that, when scaled to the DM reference frame, the 



No result of this paper is changed for Co 7^ 1- 
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Figure 7. The universal rotation curve of spiral galaxies at different luminosities and velocities [panels (a) and (b), respectively]. Radii 
are in units of R20O1 the radius encompassing a mean halo overdensity of 200, which represents the characteristic scalelength of the DM 
distribution. 



halo contribution to the RCs are essentially self-similar. Only one parameter (e.g., the total mass) specifies completely the 
halo properties (e.g., the velocity amplitude) far from the region where the LM has collapsed. 

This evidence (cf. Fig. 4) disproves the concept of 'featureless flat RC and the resulting prejudice of a 'cosmic cospiracy', 
according to which some physical process would have fine-tuned the structural parameters of the dark and visible matter, 
forcing every RC to be flat and featureless (e.g., Bahcall & Casertano 1985). Given the favour often granted to such views 
(but see Salucci & Frenk 1989 and Casertano & van Gorkom 1991), we want to stress that the large variety of RC shapes and 
the presence of clear features in some of them imply that the physical quantities underlying the structure of spirals (e.g., the 
dark-to-total mass ratio) take a wide range of values, though they are inter-related. 

We conclude this section by arguing that, since the stellar mass-to- light ratio in spirals spans a range which is limited com- 
pared with the range in luminosity, any empirical dependence on luminosity may actually reflect a more physical dependence 
on the amount of luminous matter Mlm- 



5 DISCUSSION 

• Where is DM? It is generally believed that the dark matter dominates the outer regions of spirals, while the LM dominates 
in the innermost ones. This is only partially true: in low-luminosity galaxies the transition between the two regimes occurs 
well inside the optical radius. Let us define the transition radius Rt as the innermost radius at which the LM fails to account 
for the observed RC by more than 4 times the typical observational error, i.e.: Vi U m(Rt) = 0.85 V(Rt; L). Notice that for 
R < Rt, Vium — V . So Rt is the innermost radius at which we unambiguously detect the presence of dark matter. We also 
compute the fraction Jlm of the total luminous mass inside Rt: for an exponential thin disc, /lm = 1 — (1 + x t )e~ xt with 
x t = 3.2Rt/R op t- In Fig. 8 we plot R t and Jlm as a function of luminosity. We find that in spirals of low luminosity the mass 
discrepancy begins before the disc half-mass radius xt — 1/2, so that the DM phenomenon is present also in the region where 
most of the ordinary matter resides. Conversely, at high luminosities, the phenomenon is evident only where the luminous 
matter is almost absent. 

• The Nature of DM. We are able to investigate, with a sufficiently large sample, the properties of DM in the region where it 
dominates. In Fig. 9 we show the local mass-to- light ratio dM J^ / d L JJp for galaxies of different luminosities. For our purposes 
we can assume M(R) = G~ 1 V 2 R while the luminosity distribution is obtained from the surface brightness profiles. The local 
M/L in the B-band easily reaches the value of 100, implying that the DM cannot be made of low-mass hydrogen-burning 
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Figure 8. The innermost radius Rt (in units of R op t where the DM is detected dynamically (left-hand panel) and the fraction of luminous 
matter (or, equivalently, of luminosity) inside Rt (right-hand panel) arc plotted versus luminosity. (L„ corresponds to = —20.5 and 
M| = -21.9.) 
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Figure 9. Local dynamical mass-to-light ratios (in solar units) as a function of radius at different luminosities. 



stars. Obviously, this argument does not rule out other dark baryonic candidates such as black-holes and jupiters, but obliges 
them to (i) be distributed very differently from the luminous baryonic matter, and (it) coexist with the luminous matter at 
about ~ (1/2) Ro P t- 

• The Final State of the Baryonic Infall. Flores et al. (1993) have proposed a model for the infall of pre-galactic gas into DM 
halos which does not predict, for V op t^ 100 km s _1 , any relationship between the luminosity and the RC slope, but rather 
predicts a random variation around the value of V pre d ~ O.lto.2- Flores et al. wondered whether the trend between V and 
luminosity found by PS91 may arise as a spurious effect of the small size of the sample used, of the limited extent of some RCs 
and of the assumption of an exponential thin disc taken to derive R op t- The selection criteria adopted to construct the present 
samples bypass the above-claimed problems: in particular, our data sample is larger by a factor 15 than that in PS91, and R op t 
is derived from the actual surface photometry. Notwithstanding this (or, rather, just as a result of this), the slope-luminosity 
(velocity) relationship stands out more tightly than in PS91. We have carefully investigated the crucial region 100 < V op t/(km 
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Figure 10. The URC surface 



s _1 ) < 150, which is very well sampled here: both the individual curves and the synthetic curves show a very steep slope 
(V > V pre d — 0.2), in agreement with the global trend. The latter is very strong, and emerges clearly also in samples (as 
in PS91) with many fewer (~ 50) objects and lower observational accuracy (i.e., <5V ~ 0.1). Thus, the observed RC profiles, 
unlike the Flores et al. (1993) predictions, are generally very steep, and show a marked correlation with luminosity. Navarro 
et al. (1996) show how an infall model a la Flores et al. may naturally reproduce the RC systematics highlighted in this (and 
previous) paper(s). 



6 CONCLUSION 

In this paper we have investigated the main properties of the mass structure of spirals. To do this, we have used a very large 
number (~ 1100) of galaxy RCs, to construct: (a) a sample of 131 high-quality extended RCs; and (b) a sample comprising 
616 medium-quality RCs that, co-added in 11 (x 2) synthetic curves, have thouroughly covered the whole luminosity (velocity 
amplitude) sequence of spirals. Our analysis extends out to 2 optical radii and spans approximately 6 mag. 

Both samples show that spiral RCs follow a common pattern: their amplitudes and profiles do not vary freely among 
galaxies, but depend on luminosity. At low luminosities the RCs are steep for i£$ R op t, and grow monotonically to a probably 
asymptotic value at outer radii. At high luminosities the RCs are fiat (and even decreasing) for R< Ropt, and gently fall, from 
~ Ropt outwards, to reach a probably asymptotically constant value farther out. They are very well represented by: 

1/2 

km s _1 (14) 



0.72 + 0.44 log A) + 1.6c-°-" £ 1 



R 

^^R^? = V(Ropt) LV""" ' L.) (z 2 + 0.78 2 )i« ' — 3,2 + 1.52(^)0, 

with x = R/Ropt)- The Universal Rotation Curve in eq.(14) (see Fig. 10) describes any rotation curve at any radius with a 
very small cosmic variance. In fact eq.(14) predicts rotation velocities at any (normalized) radius with a typical accuracy of 
4%. 

On the other hand, by slicing the URC to match individual observed RCs we can derive galaxy luminosities and therefore 
measure cosmic distances with a typical uncertainty of 0.3 magnitudes. The benefits of using the URC as a distance indicator 
are discussed by Hendry et al. (1996). 

A particular feature of the Universal Rotation Curve is the strong correlation between the shape and the luminosity 
(velocity) established in previous papers and confirmed here over a factor 150 in luminosity (factor of 5 in velocity) and over 
the variety of the RC profiles. This relationship sweeps a narrow locus in the profile/amplitude/luminosity space, implying 
that the great majority of spirals belong to the same kinematical family. At high luminosity, the profiles of RCs are only 
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log(a) 



Figure 11. The curve populated by spiral galaxies in the (luminosity)-(dark/visible mass ratio)-(DM core radius)-(central DM density) 
space. (The halo velocity core radius, a, is in units of the optical radius R op t-) Different luminosities are represented as grey circles along 
the curve with areas proportional to luminosity at constant Alog(L) intervals between 0.06L* and 3L» . 



mildly discrepant from the pure LM predictions, at least within R op t, so the DM content, although detectable, is modest. At 
low luminosities, on the other hand, the discrepancy between LM predictions and RC profiles is very large: the DM is the only 
relevant mass component. The mass structure in spirals is self-similar, but the dark and luminous components are coupled 
with a luminosity-dependent strength. This picture is quantified, by means of the results of Sections 3 and 4, by the following 
relationships (see Salucci et al., in preparation): 
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(where Qbbn is the cosmological density parameter of Big Bang synthesized baryons). 

A curve in the (mminosity)-(dark-to-luminous mass ratio)-(DM core radius)-(central DM density) 4-D space (see Fig. 11) 
is the geometrical analogue of these structural relationships. This curve, which is the counterpart, for spirals, of the Funda- 
mental Plane for ellipticals, shows that the distribution of DM is completely different from that of the luminous matter and 
strongly varies with luminosity, and that, especially at low luminosity, the presence of DM is already detected well inside the 
optical regions. The implications are as follows: 

• the Universe becomes darker at smaller scales; 

• smaller discs have lower densities and mass-to-light ratios; 

• smaller discs have lost a larger fraction of their original baryon content; 

• smaller haloes have higher central densities; 

• the formation of spirals is controlled by just one quantity, e.g. the total mass (or, alternatively, the initial DM overdensity) . 
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Figure 12. The luminous-mass profiles for the galaxies in Sample B, grouped by luminosity bins. For each individual object, the light 
profile L(r) oc I(r') r' dr' is normalized to its total value Loo; the radius is normalized to R op t- Grouping the light profiles by velocity 
amplitude yields a similar result. 



Whitmore, B.C., Forbes, D., & Rubin, V.C. 1988, ApJ, 333, 542 



APPENDIX A 

In this Appendix we study the connection between the velocity gradients V and 8 and the distribution of DM. We start 
from the condition of centrifugal equilibrium: 

V\R) = VZmiR) + V h 2 (R), (Al) 

where V{^ rn (R) is the quadratic sum of the three luminous components, i.e. gas, disc, and bulge. The gas contribution 
is obtained from the HI surface density, and the functional form of the stellar contribution is determined by the surface 
brightness profile (the effect of the bulge is studied in Persic, Salucci & Ashman 1993). Let us define 

, _ Vium (2R opt ) - V lum (R opt ) 

Vlum(-tiopt ) 
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Figure 13. V (um and 

^lum as a function of velocity. The data are taken from Begeman (1987), Brocils (1992a, b), Carignan & Puchc 
(1990a, b), Carignan, sancisi & van Albada (1988), Cote, Carignan & Sancisi (1991), Puchc, Carignan & Bosnia (1990), Puchc, Carignan 
& van Gorkom (1991b), and Puche, Carignan & Wainscoat (1991a). The slanting and the horizontal line represent, respectively, the 
locus populated by the observed RC slopes and the exponential disc case. 



and 



d\ogV t urn (R) 



dlogR 



(A3) 



Ropt 



and define S h and similarly for the halo component. Finally, let us define the LM fraction inside the optical radius as: 



(A4) 



V 2 (R) 

The RC-profile equation [cf. eq.(4) in Persic & Salucci 1990b] is: 

V = /3V lum + (l-/?)Vk, {AS) 

while from the definition of S we have: 

(l + S) 2 = (3(l + S lum ) 2 + (1 - (3) (1 + 5 h f . (A6) 

We investigate whether the light profile of spirals depends on luminosity by co-adding the surface brightnesses of the 616 
galaxies of Sample B to form synthetic light profiles, arranged by luminosity (velocity) bins (see Fig. 12). 

We stress that the actual light distributions do not depend on luminosity (velocity) , which implies that the velocity profile 
of the stellar component, 



V 2 (R) + V 2 (R) 
V b 2 (R opt ) + V 2 (R opt ) 



1/2 



does not depend on luminosity. In particular, at 2R opt the light has already converged: this implies that the bulge and the 
disc are already very close to a Keplerian regime, Vb(r) oc r -1 ^ 2 and Vd oc r~ 1//2 , in any object. 
We can estimate X7i um from published data taking into account that 

2 /..\2 /..\2 



Vd 

Vlurr. 



Pd + 



V b 



b + 



09 



(where the subscripts d,b,g denote the disc, bulge and gas components, respectively), and 



V 

v I'll 



Vd + V^ + Vg . 



(A7) 



(A8) 



In Fig. 13 we plot V;„ m as a function of luminosity: we find 

V lum = -0.24 ±0.03 ~ V d , (A9) 

with the residuals uncorrelated with luminosity (but weakly correlated with Hubble type). This result can be easily understood. 
In fact, (i) most of the light profiles are true exponentials for which such a result holds automatically, (ii) the discrepancies of 
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Figure 14. Synthetic rotation curves for Sample B arranged by velocity amplitude. Galactocentric radii are normalized to R op t, the 
radius encompassing 83% of the total / light. The last panel shows the rotation curve predicted for a pure self-gravitating exponential 
thin disc. 



the surface brightness profiles with respect to an exponential are independent of luminosity, (Hi) neglecting the bulge affects 
the determination of V very little (see Appendix of Persic et al. 1993), and (iv) the gas content is relatively small, and its 
effect on V; um is opposite to those considered in points (ii) and (in). 

The external gradient Si um can be estimated taking into account that 

S lum ~ -0.26 (3 d - 0.29 p b + (3 g 8 g . (AW) 

In Fig. 13 we also plot 8i um as a function of luminosity: we find that there is very little dependence of this quantity on 
luminosity, so that we can assume: 

Sium = -0.25 ±0.04. (All) 

Since V and 8 depend on luminosity while V; um and 5i U m do not, then at least two of the halo quantities 1 — /3, Vh, 8h will 
also depend on luminosity; on the other hand, a luminosity dependence of the latter quantities will generally be propagated 
into the RC profiles. 

APPENDIX B 

In this Appendix we show: 

a) the rotation curves of Sample B co-added in velocity bins (Fig. 14); 

b) the rotation curves of Sample B suitably co-added in order to probe the kinematics of the outermost regions (Fig. 15). 
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of Sample B (filled circles; the smaller symbols represent the inner curves, plotted with finer binning in Figs.l and Bl). The i 
number of luminosity /velocity bins is necessary to compensate for the sparseness of the outer data. 
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APPENDIX C 

The self-gravity of an infinitely thin disc with surface mass density distribution / = Ioe R ^ RD yields an equilibrium circular 
velocity given by: 



Vi(y) = ^/(/oifo-/A), 



(CI) 



where y — R/Rd and /„ and K„ are the modified Bessel functions computed at y/2. For mathematical simplicity, in the 
range 0.04 < x = R/R opt < 2 we approximate the r.h.s of eq.(Cl) with the much simpler function 

GM LM 1-97 x 1 - 22 



Va(R) = 1.1 



(C2) 



Rapt (x 2 + 0.78 2 ) 1 - 43 ' 

Notice in Fig. 16 that our approximation agrees with the exact function to 0.4% over the quoted radial range. For R > 2R op t 
the Keplerian regime is practically attained, and hence: 

-1/2 



V d (R) = V d (2R opt ) 



R 



2R, 



opt 



(C3) 



APPENDIX D 
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Figure 16. The approximated Vd plotted alongside the exact function (upper panel) The fractional difference is also shown (lower 
panel). 



In this appendix we present the relevant quantities of Sample A (Table 1) along with the surface-photometry and 
rotational- velocity references (Table 2). 
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N 55 
N 224 
N 247 
N 253 
N 300 
N 598 
N 628 
N 697 
N 753 
N 801 
N 891 
N 925 
N 1035 
N 1085 
N 1090 
N 1097 
N 1114 
N 1247 
N 1365 
N 1417 
N 1560 
N 1832 
N 2336 
N 2403 
N 2558 
N 2595 
N 2742 
N 2841 
N 2903 
N 2998 
N 3109 
N 3145 
N 3198 
N 3200 
N 3223 
N 3992 
N 4013 
N 4062 
N 4236 
N 4258 
N 4348 
N 4565 
N 4605 
N 4682 



Table 1. Sample A: Physical Parametei 

Dist. R„ v , M n i 

1.3 4.1 -17.83 79 

0.7 15.7 -20.81 78 

2.1 7.7 -17.52 75 
3.0 8.9 -19.82 72 

1.2 4.4 -16.83 50 



43.0 
67.7 
79.3 



24.1 
31.7 



22.2 
21.2 



12.8 
12.9 



13.7 
23.5 



-20.12 
-21.16 
-21.54 
-21.26 
-20.50 
-19.24 
-19.04 
-21.92 
-20.72 
-20.71 

-21.07 
-20.27 
-21.15 
-16.80 
-20.25 
-21.79 
-19.24 
-20.43 
-21.15 
-20.10 
-21.83 
-20.52 
-21.16 
-16.43 
-21.22 
-20.17 
-21.09 
-21.76 
-21.67 
-19.73 
-19.47 
-18.74 
-20.73 
-19.52 
-20.86 



215-039 
235-016 
240-011 



347-033 
350-023 
352-053 
374-027 
376-002 
379-006 



437-030 
439-018 
439-020 
444-047 

444- 086 

445- 058 

446- 044 
481-002 
499-005 
502-002 
507-007 
509-091 
533-004 
543-012 
548-032 
555-016 

563- 014 

564- 020 
566-022 
601-009 

M- 3- 1042 



54.3 
94.8 
37.2 
26.9 
67.4 
42.2 
65.9 



70.0 
41.4 



20.9 
20.8 
51.2 
109.2 
51.9 



16.7 
57.5 
48.2 



14.2 
14.7 
11.0 



10.5 
10.3 
12.7 
12.4 



10.8 
12.0 



10.7 
14.6 



-20.74 
-20.53 



-20.20 
-20.55 



-20.43 
-19.87 



-20.46 
-21.71 
-19.91 
-19.44 
-19.95 



-20.84 
-20.25 
-21.10 
-20.03 
-19.23 



30.0 
35.9 
21.2 



N 4800 
N 5033 
N 5055 
N 5371 
N 5585 
N 5673 
N 5905 
N 5907 
N 6503 
N 6674 
N 7083 
N 7331 
N 7339 
N 7536 
N 7591 
N 7593 
N 7606 
N 7631 
N 7793 
I 467 
I 2974 
U 2259 
U 3269 
U 3282 
U 4375 
U11810 
U12417 
U12533 
U12810 
27-008 



26.5 
22.0 



-20.67 
-20.65 
-21.T2 
-18.47 
-20.30 
-20.87 
-20.69 
-18.62 
-21.33 
-21.27 
-20.92 
-19.21 
-20.18 
-20.34 
-19.61 



-19.59 
-20.39 



-21.45 
-20.48 
-21.71 



NOTES. 

Colli}: galaxy name. 
Col\l}: morphological type. 
Col}3}: adopted distance, in Mpc. 
Col\i}: optical radius (R apt ), in kpc. 

Colfi}: absolute blue magnitude, from the Bj apparent magnitude reported in the RC3 
catalog (de Vaucouleurs et al. 1991) and the distance reported in col. (3), adopting Tully's 
(1988) correction for inclination and Burstein k Heiles's (1984) correction for Galactic 
extinction. 



Coll6}: inclinatio 



.ngle, in degrees. 



Coh.l%l&}\ RC parameter 
0.5 < R/R llT , < 1 (in km s" 



.nd Vj, such that V(R) = V + V^^) in the range 



Col.19}: velocity difference V(2R opi ) - V(R aF 



-21.04 
-20.89 
-17.75 



-19.62 
-21.01 
-21.40 
-21.40 



© 0000 RAS, MNRAS 000, 000-000 



22 



M.Persic, P.Salucci and F.Stel 



Object" 


S.Ph. 6 


RC" 


Dist. 


N 55 


lo 


16 




224 


4 


34 




N 247 




18 




253 




19 




N 300 


18 


20 


1 


N o98 


19 


11,21 




N 628 




27 




N 697 


30 


30 


10 


N 753 


1 


1 


2 


N 801 


1,28 


1 


2 


N 891 


20 


22 


8 


N 925 


10 


13,33 


8 


N 1035 


1 


1 


8 




1 


2 


3 


XT i rnn 
JN 1090 


30 


30 


10 


N 1097 


30 


5 


1 


N 1114 


30 


30 


10 


N 1247 


30 


30 


10 


N 1365 


12 


23 


8 


N 1417 


1 


2,30 


3 


1560 


29 


31 


8 


N 1832 


30 


30 


10 


N 2336 


8 


4 


8 


N 2403 


10 


6 


8 


N 2558 


26 


29 


6 


N 2595 


27 


29 


6 




1 


1 


1 


N 2841 


4 


6 


8 


XT 1 (1 (1 9 

JN 2903 


10 


6 


8 


XT 1 (1 (1 

.n 2998 


1,28 


1,32 


2 


N 3109 


21 


24 


1 


3145 


30 


2 


3 


N 3198 


4 


6 


8 


xt tiirn 
JN 3200 


1 


2 


3 


N 3223 


30 


2 


1 


N 3992 


2 


8 


8 


N 4013 


23,31 


17 


8 


N 4062 


1 


1 


8 


N 4236 


4 


3 


8 


N 4258 


4 


9 


8 


N 4348 


30 


30 


10 


N 4565 


9,24 


26 


8 


N 4605 


1 


1 




N 4682 


1 


1 


2 



Table 2. Sample A: References 



Object" 


S.Ph. 6 


RC 


Dist. 


\T A Qflfl 

JN 4oUU 


1 


2 


1 


N 5033 


10 


6 


8 


N 5055 


10 


10 


8 


N 5371 


10 


6 


1 


N 5585 


22 


25 


8 


N 5673 




7 


4 


N 5905 




12 


5 


N 5907 


23 


26 


8 


N 6503 


28 


35 


1 


N 6674 


7 


32 


9 


7083 




2 


1 




4 




8 




30 


30 


10 


N 7536 


27 


29 


6 




30 


14 




N 7593 


27 


29 


6 


JN ( bUb 


1 


2 


1 






29 




N 7793 


2o 


28 


1 


1467 


1 


1 


1 




30 


30 


10 


U 2259 


13 


15 


1 


U ozby 


27 


29 


7 




27 


29 


7 


U 4375 


14 


14 


6 


U11810 


1 


2 


3 




30 


14 


6 


T T 1 1 T. 9 9 

U 12033 


30 


30 


10 


TT 1 10 1 (1 

U 12810 


1 


2 


3 


TT IT AAO 

27-008 


30 


30 


10 


TT 9n irnn 
lj 30-009 


30 


30 


10 


TT a ft nn 
lj 40-012 


30 


30 


10 


TT 11 nnn 
lj 41-009 


30 


30 


10 


tt a a n i i 
li 09-011 


30 


30 


10 


E 71-005 


30 


30 


10 


E 75-037 


30 


30 


10 


E 82-008 


30 


30 


10 


E 88-016 


30 


30 


10 


E116-012 


30 


30 


10 


E121-006 


30 


30 


10 


E123-023 


30 


30 


10 


E141-020 


30 


30 


10 


E141-034 


30 


30 


10 


E184-051 


30 


30 


10 



Object" 


S.Ph. 6 


RC" 


Dist. 


TM 1 C 11 9(1 

fi21O-039 


30 


36 


10 


TT i 9 t. n i a 
fi230-01b 


30 


30 


10 


Tallin (111 
fi240-011 


30 


30 


10 


TT o a a n i a 
fi2b9-019 


30 


30 


10 


E282-003 


30 


30 


10 


TTIO.1 (11/1 

fi284-024 


30 


36 


10 


TT o o a nic 

li28b-01b 




30 


10 


TT 1 T (119 

fi28 (-013 


30 


30 


10 


E289-010 


30 


30 


10 


TTH1I1 (1(1/1 

1i299-004 


30 


30 


10 


TT 9 (Ifi (191 

li30b-032 


30 


30 


10 


TT9T1 (1/1 C 

fi322-040 


30 


36 


10 


^1911 117(J 

rj322-07b 


11 


36 


10 


E346-014 


30 


30 


10 


TT 9 1 T (19 9 

fi34 (-033 


30 


30 


10 


TT9C(1 (119 

fi3o0-023 


30 


30 


10 


TT9 CI (1 C9 

ri3o2-0o3 


30 


30 


10 


E374-027 


30 


30 


10 


TT9TC (1(11 

iiA (b-002 


30 


f 


10 


TT9TO (l(lfl 

li'3 /9-00b 


30 


30 


10 


TT 9 9 (1(11 

ii383-002 


30 


36 


10 


E383-088 


30 


36 


10 


TT 1 9 7 (1 9(1 

i/43 (-030 


30 


36 


10 


TTM 911 (110 

Ii439-018 


30 


36 


10 


TT 1 9 Ci (11(1 

f/439-020 


30 


36 


10 


ri 1 1 d Hi? 
E444-047 


30 


36 


10 


TP A A A (1 0(1 


30 


36 


10 


TTM A C (ICO 

ri44o-0o8 


30 


36 


10 


TT 1 1 flA A 

i;44b-044 


30 


30 


10 


El 1 O 1 (1(11 

E481-002 


30 


30 


10 


TT/IIKl (1(1C 

ii499-000 


30 


30 


10 


TT I (11 (1(11 

iiO02-002 


30 


36 


10 


TT T. (IT (KIT 

iiOO (-00 ( 


30 


30 


10 


E509-051 


30 


36 


10 


E533-004 


30 


30 


10 


E543-012 


30 


30 


10 


E548-032 


30 


30 


10 


E555-016 


30 


30 


10 


E563-014 


30 


30 


10 


E564-020 


30 


30 


10 


E566-022 


30 


30 


10 


E601-009 


30 


30 


10 


M-3-1042 


30 


30 


10 



Galaxy name. 

Surface photometry references: 
1) Kent 1986; 2) Elmegreen & Elmegieen 1985; 3) Shostak & van der Kruit 1984; 4) Kent 1987; 5) Blackman 
& van Moorsel 1984; 6) Kent 1985; 7) Broeils & Knapen 1991; 8) Boroson 1981; 9) van der Kruit & Searle 
1981a; 10) Wevers, van der Kruit & Allen 1986; 11) Schommer et al. 1993; 12) Hodge 1978; 13) Carignan, 
Sancisi & van Albada 1988; 14) Whitmore, Forbes & Rubin 1988; 15) Puche, Carignan & Wainscoat 1991; 

16) Carignan & Puche 1990b; 17) Puche, Carignan & van Gorkom 1991; 18) Puche, Carignan & Bosma 1990; 
19) Regan & Vogel 1994; 20) van der Kruit & Searle 1981b; 21) Jobin & Carignan 1990; 22) Cote, Carignan 
& Sancisi 1991; 23) van der Kruit & Searle 1982; 24) Shaw & Gilmore 1989; 25) Carignan & Puche 1990a; 
26) Schombert & Bothun 1987; 27) Amram 1992, private communication; 28) Bottema 1989; 29) Broeils 
1992a; 30) Mathewson, Ford & Buchhorn 1992; 31) Bottema, Shostak & van der Kruit 1987; *) No surface 
photometry available: R op t = with #25 taken from RC3. 

Rotation velocity references: 
1) Rubin, Ford & Thonnard 1980; 2) Rubin et al. 1982; 3) Shostak 1973; 4) van Moorsel 1983b; 5) Ondrechen, 
van der Hulst & Hummel 1989; 6) Begeman 1987; 7) van Moorsel 1983a; 8) Gottesman et al. 1984; 9) van 
Albada 1980; 10) Bosma 1981a; 11) Newton 1980; 12) van Moorsel 1982; 13) Gottesman 1980; 14) Rubin, 
Whitmore & Ford 1988; 15) Carignan, Sancisi & van Albada 1988; 16) Puche, Carignan & Wainscoat 1991; 

17) Bottema, Shostak & van der Kruit 1987; 18) Carignan & Puche 1990b; 19) Puche, Carignan & van 
Gorkom 1991; 20) Puche, Carignan & Bosma 1990; 21) Rogstad, Wright fc Lockhart 1976; 22) Sancisi & 
Allen 1979; 23) Ondrechen & van der Hulst 1989; 24) Jobin & Carignan 1990; 25) Cote, Carignan & Sancisi 
1991; 26) Sancisi & van Albada 1987; 27) Kamphuis & Briggs 1992; 28) Carignan & Puche 1990a; 29) 
Amram et al. 1992; 30) Persic & Salucci 1995; 31) Broeils 1992a (approaching half only); 32) Broeils 1992b; 
33) Marcelin, Boulesteix & Courtes 1982; 34) Kent 1989; 35) Bottema 1989; 36) Schommer et al. 1993. 

Distance references: 

1) Tully 1988; 2) Rubin, Ford & Thonnard 1980; 3) Rubin et al. 1982; 4) van Moorsel 1983a; 5) van Moorsel 
1982; 6) Rubin, Whitmore & Ford 1988; 7) Amram et al. 1992; 8) Burstein 1991 (MARK II release); 9) 
Broeils 1992b; 10) Recession velocity measured in Galactic Standard of Rest (Vcsr), from RC3; *) estimated 
from the ff-band Tully-Fisher relation (Aaronson et al. 1982), with the ^-magnitude from Whitmore 1984. 
Distances quoted from refs. 2) through 6) have been scaled to a Hubble constant value of 75 km s -1 Mpc -1 , 
and distances computed from refs. 8) and 10) likewise assume H a — 75 km s _1 Mpc -1 . 
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